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ABSTRACT 


Continuing  wavelength  dependent  absorption  measure¬ 
ments  of  "high  purity"  GaAs  grown  from  the  melt  by  a  variety 
of  techniques  show  remarkably  similar  structure  in  all  samples. 
There  appears  to  be  a  definite  lower  limit  of  about  3“  .008  cm 
at  10.6  pm  for  melt  grown  GaAs.  Wide  band  emissivity  measure¬ 
ments  performed  by  Stierwalt  (NELC  San  Diego)  agree  well  with 
our  data,  except  for  a  small  shift  of  scale  attributable  to  a  dif¬ 
ference  in  ambient  temperature.  Preliminary  measurements  of 
the  dielectric  constant  of  GaAs  have  been  obtained.  A  systematic 
study  of  diffuse  intensity  streaks  in  the  electron  diffraction  pattern 
of  GaAs  single  crystals  indicates  that  they  are  of  thermal  origin, 
and  do  not  arise  from  inclusions  or  impurities.  The  absorption 
coefficient  of  a  good  quality  sample  of  GaAs  has  been  reduced  to 
that  of  best  quality  samples  by  annealing.  Greatest  reduction 
occurred  in  the  9.  5|im  feature  common  to  all  melt  grown  samples 
measured.  High  purity  KC1  crystals  have  been  grown,  and 
preparations  are  now  being  made  for  the  growth  of  OH  doped  KC1 
and  high  purity  KBr.  Experiments  on  self  diffusion  of  Cd  and  Te, 
and  high  temperature  Hall  effect  in  CdTe  and  In-doped  CdTe  have 
been  completed.  Experiments  have  been  performed  to  determine 
10.6  Um  absorption  in  a  highly  absorbing  material  (quartz)  using 
acoustical  surface  wave  techniques.  The  data  require  numerical 
analysis  (now  in  progress)  to  extract  the  absorption  coefficient. 

An  analysis  of  the  influence  of  stress  induced  birefringence  on 
thermal  slewing  of  acentric  beams  has  been  substantially  completed. 
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1. 


INTRODUCTION 


The  format  of  this  report  follows  closely  that  of  the  first  quarterly 
report  in  which  projects  are  identified  by  codes  keyed  to  the  contract 
work  statement. 

The  various  categories  are  briefly: 

a)  Crystal  growth 

b)  Polycrystalline  window  fabrication 

c)  Mechanical  properties  of  window  materials 

d)  Window  material  defect  characterization 

e)  Theory  of  residual  IR  optical  absorption 

f)  Absorption  measurement  techniques 

g)  Theoretical  evaluation  of  optical  performances  of 
windows 
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a.  1  Effect  of  Oxygen  and  Other  Impurities  on  IR  Absorption  in  II- VI  and 
III- V  Compounds 

James  M.  Whelan 

A.  GaAs  Epitaxial  Growth- Horizonal  Reactor 

The  apparatus  has  been  completed  and  is  presently  being  outgassed. 

It  will  be  available  during  the  next  quarter.  Screening  examinations  were  in 
part  responsible  for  the  fact  that  film  growth  was  delayed. 

B.  GaAs -Oxygen  Studies 

Progress  was  delayed  because  of  a  (successful)  Ph.  D.  qualifying  exami¬ 
nation.  Most  of  the  work  this  quarter  was  spent  on  evaluating  the  zirconia 
electrolytic  cell  as  a  means  of  establishing  the  relationship  between  the  fugacity 
of  oxygen  as  a  function  of  the  atomic  fraction  of  oxygen  dissolved  in  gallium  at 
500°C.  This  information  is  needed  to  develop  a  means  for  measuring  oxygen 
in  GaAs  and  other  materials.  Difficulties  in  identifying  the  extent  of  hydrogen 
ion  diffusion  in  the  zirconia  have  introduced  an  uncertainty  in  the  above  relation¬ 
ship.  Our  preliminary  one  at  500°C  is: 

—  -11  -  — 

P2(Oz)  -  1.4x10“  X(O)  atm  2  where 

P(0,)  is  the  effective  oxygen  pressure  in  equilibrium  with  a  gallium  solution 

&  -6 
in  which  the  atomic  fraction  of  oxygen  is  X(O).  This  is  f  r  values  of  X(O)  s10 

Preliminary  indications  are  that  the  value  of  the  constant  is  not  strongly  tempe¬ 
rature  dependent  and  that  uncertainties  in  the  interpretation  of  cell  errif  data  are 
substantially  less  at  600°C. 


a.  2  Optimization  of  Alkali  Halide  Window  Materials 


Paul  J .  Sulicfata,  R.  E.  Chaney 

Single  crystals  of  KC1,  purified  by  ion- exchange  ai  d  treatment  with  HC1 
and  C^.  have  been  grown  by  the  procedures  described  in  previous  reports.  The 
infrared  absorption  of  these  crystals  is  now  being  studied  by  laser  calorimetry 
and  low-temperature  spectrophotometry;  results  will  be  tabulated  in  the  next 
report. 

Preparations  are  now  underway  for  the  growth  of  KC1  crystals  containing 
OH  ions.  Two  methods  of  doping  will  be  used;  (1)  introducing  KOH  into  the 
molten  KC1  and  (2)  exposing  the  molten  salt  and/or  grown  crystals  to  water 
vapor.  These  cyrstals  will  also  be  studied  by  calorimetry  and  infrared  spectro¬ 
photometry. 

Preparations  are  also  being  made  for  the  growth  of  purified  KBr  crystals. 
Additional  ion-exchange  columns  have  been  assembled  and  gas-handling  systems 
have  been  installed  for  treatment  of  the  molten  KBr  with  HBr  and  Br,,.  This 

w 

facility  should  be  operational  sometime  during  the  next  quarter. 
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GROWTH  OF  CRYSTALS  FOR  XR  WINDOW  RESEARCH 


a.  3 


E.  A.  Miller  and  W.  R.  Wilcox 
(Chemical  Vapor  Deposition  of  Bulk  Gallium  Arsenide) 

Work  during  this  quarter  was  directed  toward  growth  of  bulk  gallium 
arsenide  by  the  open  tube  chemical  vapor  deposition  method  using  AsCl^, 
gallium  and  hydrogen.  Results  are  summarized  in  Table  I. 

o 

The  all-glass  system  with  furnace  temperature  control  to  +  2  C  was 

placed  in  operation  on  schedule.  Substrates  used  in  this  series  of 

experiments  were  GaAs  oriented  (100)  2°  off  toward  (110)  and  were 

17  18 

tellurium  doped  in  the  range  of  5  x  10  to  2  x  10  atoms.  The  usual 
cleaning  and  preparation  procedure  was  as  follows:  The  wafers  were 
mechanically  polished  to  a  bright  mirror  finish  and  then  lightly  etched  in 
5  si  si  H2S0^sH202:H20  for  10  minutes  to  remove  traces  of  polishing  com¬ 
pounds  and  damage.  Then  substrates  were  rinsed  in  H20,  methanol,  and 
dried  in  a  stream  of  nitrogen.  The  substrates  were  then  etched  at  900°C 
in  the  reactor  just  before  deposition  at  740-780°C. 

The  first  few  runs  accomplished  only  etching  of  the  substrate. 
Attempts  were  made  to  modify  the  vapor  composition  in  the  region  of  the 
substrate  in  the  hope  of  arriving  at  the  proper  conditions  for  growth. 

One  possible  solution  involved  a  modified  Ga  source  holder  consisting  of  the 
thermocouple-well  extending  only  to  the  upstream  end  of  the  Ga  melt.  This 
should  have  allowed  much  better  vapor  contact  with  the  Ga  melt.  However, 
negative  results  were  obtained.  Another  modification  was  placer  ent  of 
quartz  discs  in  the  tube  near  the  Ga  source  and  substrate  to  cut  down  on 
free  convection  of  the  depositing  species.  This  resulted  in  an  amporhous 
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TABLE  1.  SUMMARY  OF  GaAs  CVD  RUNS 
(3/1/73  -  5/31/73) 
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appearing  GaAs  deposit.  Doubling  the  gas  flow  rate  had  no  noticeable 
effect. 

Deposition  conditions  on  {lOO}  oriented  substrates  proved  to  be 

much  more  critical  than  on  previously- tested  {111}  A  and  B  substrates. 

Therefore  we  proceeded  to  modify  the  system  by  obtaining  a  longer  furnace 

and  establishing  a  temperature  gradient  that  was  reported  to  give  good 

* 

results  on  {100}  substrates  with  this  AsCl^#  Ga  system.  Figure  1  shows 
the  temperature  profile  established  in  this  furnace.  This  new  system  had 
a  smaller  diameter  tube  (40  mm  Z.D.)  and  an  open  boat  inserted  into  the 
tube  for  the  gallium  source.  The  hydrogen  flow  through  the  AsCl3  bubbler 
was  usually  120  cc/min  and  the  hydrogen  dilution  flow  was  also  120  cc/min. 

A  palladium-diffusion  hydrogen-purification  system  was  also  installed.  A 
temperature  controller  held  the  temperature  of  the  AsC13  bubbler  to  +  0.1°C. 
In  subsequent  runs  we  located  the  substrate  in  the  tube  near  the  best 
appearing  GaAs  wall  deposits  from  preceding  runs,  although  this  position 
shifted  slightly  from  run  to  run.  Finally,  we  were  successful  and  obtained 
about  ?0  urn  of  growth  on  a  {100}  substrate  in  24  hours.  Many  growth 
hillocks  were  found  on  the  surface. 

We  achieved  better  results  when  the  substrates  were  placed  under  a 
cover  during  saturation  of  the  gallium  melt  prior  to  deposition.  We  be¬ 
lieve  that  without  a  cover  GaCl  disproportionates  on  the  substrate  to 
GaCl^  and  metallic  gallium.  The  Ga  on  the  substrate  then  acts  as  a  pre¬ 
ferred  site  for  hillock  and  whisker  formation  by  a  vapor-liquid-solid 

* 

Wolfe,  C.M.,  and  Foyt,  A.G.,  Electrochem.  Tech.,  6,  208  (1968). 
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FURNACE  TEMPERATURE 


DISTANCE  ALONG  FURNACE  (inches) 


Figure  1.  Temperature  profile  in  GaAs  CVD  system 
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growth  mechanism.  In  Figure  2  a  portion  of  the  substrate  in  a  hillock 
free  region  is  shown.  Figure  3  and  4  demonstrate  typical  hillock  for¬ 
mation  on  these  { 100}  substrates.  (Note  the  four-fold  symmetry.) 

These  surfaces,  in  spite  of  the  imperfections,  appear  highly  specular. 

The  next  24-hour  run  consisted  of  a  {ill}  A  and  B  twin  substrate 
along  with  a  {l00}  substrate.  The  growth  rate  on  the  {ill}  A  surface 
was  much  larger  and  rougher  than  on  the  {ill}  B  and  {100},  as  expected. 
Figure  5  shows  the  result  of  this  run.  Note  the  uniformly  rough  but 
hillock-free  surface,  approximately  250  pm  thick.  This  growth,  even 
though  rough,  would  be  suitable  for  absorption  measurements  since  such  a 
rough  surface  could  be  easily  polished  mechanically. 

Next  we  attempted  to  produce  a  deposit  1  to  2  mm  thick  over  five 
days  on  a  {ill}  A  substrate.  The  results  were  poor  -  slight  etching  of 
the  substrate  was  even  observed.  An  overgrowth  of  large  dendrites 
covered  the  area  above  the  substrate.  We  also  noted  large  periodic 
fluctuations  (~30°C)  of  substrate  temperature  every  12  to  13  hours, 
probably  due  to  line  voltage  fluctuations.  The  wall  deposits  appeared 
better  upstream  from  the  substrate.  This  has  been  observed  before  —  when 
the  substrate  is  placed  in  the  region  previously  showing  good  wall 
deposits,  the  wall  deposits  now  appear  better  in  another  region. 

We  plan  to  install  constant  voltage  transformers  for  the  power  con¬ 
trollers  to  eliminate  the  temperature  fluctuation  problem.  A  longer 
substrate  holder  and  cover  that  will  hold  up  to  five  samples  is  being 
built.  Hopefully  under  any  given  set  of  conditions  a  deposit  suitable  for 
absorption  measurement  can  be  grown  thereby. 

10 


Figure  3.  Hillocks  on  deposit  on  (100)  suhstrate. 
CVD-15.  50  X. 


Figure  4.  Scanning  electron  micrograph  of  CVD-16. 
200  X. 


Figure  5.  Deposit  on  {ill}  A  and  B  twinned  substrate, 
showing  greater  deposit  thickness  on  {ill} 

A  face  on  left. 
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b.  1.  Fabrication  of  Polycrystalline  IR  Window  Materials 
S.  M.  Copley,  J.  M.  Whelan  and  V.  Rana 
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Hot  Press 

The  hot  press  and  its  ope  ration  have  been  described  in  the  previous 
quarterly  reports.  The  hot  press  has  been  cor.structed  and  tests  are  underway 
to  prepare  it  for  operation. 

As  explained  in  earlier  reports,  the  use  of  liquid  arsenic  as  a  sintering 
aid  in  the  hot  pressing  of  GaAs  powder,  is  made  possible  by  containing  it  in 
the  die  with  a  glass  seal  backed  up  by  a  suitable  pressure  of  argon  gas.  The 
use  of  an  inert  atmosphere  is  required  to  prevent  the  graphite  die  and  the 
molybdenum  chill  from  oxidizing.  In  order  to  evacuate  the  hot  press  chamber 
before  introducing  the  argon  gas,  the  chamber  is  connected  to  a  vacuum  system 
consisting  of  a  mechanical  pump  and  a  cold  trap.  It  is  inevitable  that  there  are 
a  large  number  of  joints  present  ir.  the  system.  After  a  considerable  amount  of 
leak  testing  and  repairing,  it  is  now  possible  to  maintain  a  vacuum  of  70  gi  of 
Hg  or  better  in  the  hot  press  chamber  with  the  mechanical  pump  isolated. 
Considerable  time  has  also  been  spent  this  period  calibrating  thermocouples 
and  determining  temperature  profiles  in  the  die  obtained  by  adjustment  of  the 
flow  rate  of  cooling  water  to  the  molybdenum  chill. 

Since  this  hot  press  can  be  used  for  compacting  powders  of  semiconductor 
materials  other  than  GaAs,  initial  experiments  have  been  planned  for  compacting 

V 

CdTe  powder.  CdTe  has  been  obtained  for  this  purpose  from  II-VI,  Incorporated. 
In  this  case  liquid  Cd  will  be  utilized  as  a  sintering  aid.  Unlike  GaAs,  compaction 
of  CdTe  can  be  done  under  one  atmosphere  pressure  of  the  argon  gas. 

Powders 

One  of  the  major  difficulties  in  producing  a  transparent  compact,  is  obtaining 
powders  of  the  semiconductor  mateiials  in  the  submicron  size  range.  In  the 
initial  experiments,  which  are  planned  only  to  check  the  efficiency  of  the  sintering 
process,  ball  milled  powder  will  be  used. 

In  experiments  carried  out  so  far,  we  have  observed  a  considerable  ?.mount 
of  iron  introduced  into  the  powder  while  ball  milling  in  a  steel  jar.  In  the  case 
of  GaAs  this  iron  can  be  leached  out,  as  GaAs  dissolves  relatively  slowly  in 
HC1.  This  operation  can  not,  however,  be  performed  in  the  case  of  CdTe  which 
* - 

II-VI  Incorporated,  Glenshaw,  Pa.  ,, 
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dissolves  rapidly  in  HC1. 

There  is  evidence  that  the  steel  introduced  into  the  powder  comes 
mostly  from  wearing  of  the  jar  rather  than  the  balls.  Hence  it  is  possible 
to  avoid  this  contamination  by  using  a  polyethylene  bottle  as  a  jar  while 
milling  the  CdTe  crystals. 

While  admittedly,  current  procedures  are  not  conducive  to  high  purity 
in  the  case  of  CdTe,  a  high  purity  process  of  general  applicability  to  III-V 
and  II- VI  compounds  is  being  explored.  In  this  method,  small  crystals  (  of 
submicron  size  range  )  are  grown  by  condensation  from  vapor  phase,  in 
an  evaporator  which  is  filled  to  a  pressure  of  a  few  tenths  of  a  torr  with  purified 
inert  gas.  With  this  residual  pressure  in  the  system,  atoms  of  semiconductor 
in  the  vapc  r  phase  have  a  small  mean  free  path.  They  collide  with  inert  gas 
atoms,  thereby  losing  energy  and  condensing  in  the  cold  parts  of  the  evaporator. 
These  condensed  atoms  will  serve  as  nuclei  for  the  growth  of  small  crystals. 
Any  resulting  agglomerate  with  loosely  adhering  particles  can  be  readily 
broken  up  by  ultrasonic  vibration. 
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c.  1.  Mechanical  Behavior  of  III-V  and  II- VI  Compounds 


S.  M.  Copley  and  V.  Swaminathan 

The  objective  of  this  investigation  is  an  understanding  of  the  mechanical 
behavior  of  various  candidate  window  materials.  Such  an  understanding  is 
needed  because: 

(1)  The  window  is  a  load-bearing  structure  which  is  subjected  to  static 
and  aerodynamic  loads. 

(2)  Thermomechanical  techniques  are  a  promising  approach  for  window 
fabrication. 

In  previous  reports  we  have  described  the  results  of  stress-strain  experiments 
on  single  crystals  of  Si-doped  GaAs.  This  work  is  continuing  and  £.n  abstract 
describing  it  has  been  submitted  to  the  American  Ceramic  Society  for  presen¬ 
tation  at  the  Fall  Meeting  at  Pittsburgh,  Pa.  .  In  this  report,  we  shall  describe 
the  results  of  experiments  to  determine  the  effect  of  various  thermal  and 
thermomechanical  treatments  on  the  microstructure  and  absorption  coefficient 
of  GaAs.  We  have  also  investigated  the  microstructure  of  several  types  of 
GaAs  with  low  absorption  coefficients  in  the  as -received  condition  and  the 
results  of  this  investigation  will  also  be  included. 

Characterization  of  As-Received  Material 

Transmission  electron  photomicrographs  have  been  obtained  for  the 
high  resistivity,  undoped  GaAs  (WA1000)  and  the  Cr-doped  GaAs  (  Bell  and 
Howell  ),  both  in  the  as-received  condition.  The  undoped  material,  Fig.  la, 
was  observed  to  have  a  very  low  dislocation  density  (  about  one  dislocation 
per  foil  ).  The  Cr-doped  material.  Fig.  lb,  had  a  somewhat  higher  dislocation 
density.  No  precipitates  were  observed  in  foils  of  either  material.  The 
dark  regions  in  Fig.  lb  are  artifacts  at  the  foil  surface 
Thermally  Diffuse  Scattering 

Diffuse  intensity  streaks  have  been  observed  in  transmission  electron 

* 

What  appeared  to  be  precipitates  were  observed  in  the  first  foils  of  the  Cr-doped 
material  that  we  examined.  We  have  been  unable  to  reproduce  these  observations 
after  examining  many  additional  foils  prepared  from  the  Cr-doped  material. 

Work  described  in  this  section  was  carried  out  in  collaboration  with  Prof.  G.H. 
Narayanan,  Dept,  of  Materials  Science,  University  of  Southern  California. 
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Fig.  2  Diffuse  intensity  streaks  in  electron 

diffraction  pattern  of  GaAs  (  electron 
beam  approximately  parallel  to  i_001]  ) 
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diffraction  patterns  of  both  undoped  and  doped  GaAs.  These  streaks,  which  are 
shown  in  Fig.  2,  can  be  regarded  as  arising  due  to  the  intersection  of  { 1 1 0 } 
sheets  of  intensity  in  reciprocal  space  with  the  Ewald  sphere.  A  systematic 
study  of  the  variation  of  streak  intensity  with  the  crystallographic  orientation 
of  the  incident  beam  indicates  that  these  streaks  can  not  be  caused  by  planar 
defects  such  as  thin  platelike  precipitates,  twins  or  stacking  faults.  Futhermore, 
no  indication  of  rod-like  precipitates  has  been  found  either  in  bright  field  or 
by  dark  field  imaging  of  the  diffuse  streaks.  It  is  believed  that  the  diffuse 
intensity  streaks  are  of  thermal  origin  similar  to  diffuse  streaks  observed  in  a 
variety  of  other  materials  including  Al,  Si,  Ge,  Fe  and  KCl  (  see  review  paper 
by  Honjo  et  al*  ).  We  have  recently  found  that  by  hot-stage  electron  microscopy 
the  intensity  of  the  diffuse  streaks  increases  with  increasing  temperature  thus 
supporting  the  hypothesis  that  the  streaks  are  of  thermal  origin.  A  paper 
describing  these  observations  is  in  preparation. 

Annealing  Treatment 

A  comparison  of  the  absorption  coefficient  of  the  Cr-doped  GaAs  (  Bell 

and  Howell  )  to  that  of  the  high- resistivity  GaAs  (WA-1000),  see  Section  f.  1 

Fig.  1  ,  suggests  that  thw.  difference  may  be  due  to  absorption  by  a  small 

2 

volume  fraction  of  precipitates.  According  to  Sparks  and  Duthler  ,  the 
absorption  cross  section  for  a  spherical  metallic  inclusion  is  given  by  the 


equation 


where 


_  ,,  J  ,  iju  .  Ta  2 

0  ,  =  12  *_  _  (  — )  —  Tra 

abs  H  uu  c 

P 


€  =  dielectric  constant  of  matrix 

H 

uu  =  frequency  of  radiation 

uu  =  plasma  frequency 

P 

T  =  relaxation  frequency 

c  =  velocity  of  light 

a  =  radius  of  the  inclusion 

The  contribution  of  the  inclusions  to  the  absorption  is 

A  8  =  0  ,  Nt 
abs  I 

where  Nj  is  the  number  of  inclusions  per  unit  volume.  Setting  =  12, 
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the  intensity  of  the  diffuse  streaks  increases  with  increasing  temperature  thus 
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and  Howell  )  to  that  of  the  high- resistivity  GaAs  (WA-1000),  see  Section  f.  1 

Fig.  I  ,  suggests  that  the  difference  may  be  due  to  absorption  by  a  small 

2 

volume  .  action  of  pr?fi ouates.  According  to  Sparks  and  Duthler  ,  the 
absorption  cross  section  for  a  spherical  metallic  inclusion  is  given  by  the 
equation  ^ 

0  ,  -  1 2  (  — )  —  TTa2  ( 1 ) 

abs  H  uu  c  '  ' 

P 


whe  re 
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=  dielectric  constant  of  matrix 
uu  =  frequency  of  radiation 

uUp  =  plasma  frequency 

F  =  relaxation  frequency 

c  =  velocity  of  light 

a  =  radius  of  the  inclusion 

The  contribution  of  the  inclusions  to  the  absorption  is 
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U)  =  5  x  10^  sec  * ,  T=  10 sec  *  and  a  =  100^,  we  obtain  a  ,  =7.6 

P  a.¥  .  1 

for  radiation  with  a  wavelength  of  10.  6  4.  Setting  £  =  1. 2  x  10  cm  .  see 

12  -3 

Section  f.  I  ,  Fig.  1  ,  we  obtain  N  =  1.58  x  10  cm  .  Since  the  volume 

1  -13 

of  material  shown  in  a  typical  photomicrograph  is  1  4  x  1  4  x  0.  1  4  =  10  cm, 

the  number  of  particles  contained  in  such  a  volume  for  the  observed  6  6  is 
0.  158  (  i.e.  ,  one  inclusion  should  be  observed  in  seven  such  volumes  ).  If  the 


absorption  is  caused  by  larger  inclusions,  they  should  have  a  lower  density. 

Thus,  if  the  absorption  difference  between  the  Cr-doped  and  undoped  materials 
at  10.  6  4  is  d\ie  to  inclusions,  these  inclusions  might  have  escaped  detection 
due  to  our  limited  number  of  observations.  It  should  be  noted  that  the  increase  in 


the  absorption  difference  between  the  Cr-doped  and  undoped  materials  with 
decreasing  wavelength  is  consistent  with  inclusion  absorption  according  to 
Eqns  (1)  and  (2). 

To  further  investigate  the  possibility  that  inclusions  were  responsible 


for  the  difference  in  absorption  between  the  Cr-doped  and  the  undoped  materials, 
it  was  decided  to  anneal  and  quench  the  Cr-doped  material.  If  the  difference  in 
absorption  was  due  to  nclusions  it  was  reasoned  that  it  might  be  reduced  by  a 
solutionizing  anneal.  Accordingly,  the  Cr-doped  material  was  sealed  in  a 
silicon  tube,  annealed  at  1100°C  for  15  minutes  and  then  rapidly  cooled.  A 
comparison  of  the  absorption  coefficients  of  the  annealed  and  as-received 
specimens  is  shown  in  Fig.  3.  The  annealing  treatment  produces  a  decrease  in 
absorption  coefficient  and  the  amount  of  the  decrease  becomes  greater  with 
decreasing  wavelength.  The  absorption  coefficient  of  the  annealed  Cr-doped 
GaAs  is  comparable  to  the  as-received  high  resistivity  undoped  (WA  1000)  material. 


Figure  4  shows  an  electron  photomicrograph  of  the  annealed  Cr-doped 
material.  The  most  obvious  change  upon  annealing  (  see  Fig.  lb  )  is  the 
reduction  in  dislocation  density. 


Cold  Working  and  Process  Annealing 

Two  approaches  may  be  distinguished  for  fabricating  windows  by  thermo¬ 
mechanical  means.  In  the  first,  called  "hot-working",  the  material  is  deformed 


at  a  temperature  that  is  sufficiently  high  for  recrystallization  to  occur.  This 
method  provides  an  attractive  approach  for  fabricating  windows  of  the  alkali 
halides  due  to  their  strong  tendancy  to  cleave  and  their  lack  of  five  independent 
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slip  systems.  In  the  second,  called  "  cold-working  and  process  annealing  ", 
the  material  is  deformed  at  a  temperature  too  low  for  recrystallization  to  occur. 
The  material  is  then  annealed  at  a  temperature  sufficient  to  effect  recrystallization. 
This  approach  is  attractive  for  GaAs  because  process  annealing  can  be  carried 
out  in  a  controlled  atmosphere  to  suppress  oxidation  and  arsenic  loss.  Further¬ 
more,  GaAs  does  not  have  as  strong  a  tendancy  to  cleave  as  the  alkali  halides 
and  the  primary  a/2  {lll}<liO>  slip  systems  in  this  material  c omprise  five 
independent  slip  systems.  Thus  substantial  ductility  at  temperatures  below 
the  recrystallization  range  is  anticipated. 

Experiments  to  explore  the  possibility  of  fabricating  GaAs  windows 
by  cold  working  and  process  annealing  have  been  carried  out  on  Si-doped  GaAs. 
Although  this  material  is  unsuitable  for  windows,  it  is  believed  that  its  response 
to  thermomechanical  treatments  should  be  the  same  as  materials  that  are 
suitable  such  as  the  Cr-doped  or  the  high  resistivity  undoped  GaAs. 

Figure  5  shows  the  dislocation  substructure  of  a  Si-doped  GaAs  crystal 
deformed  3  pet  at  550°C.  One  particularly  noticable  feature  is  the  occurance  of 
many  long  straight  segments  of  screw  dislocation.  Considerably  greater  strain 
than  3  pet  is  required  to  produce  recrystallization  upon  subsequent  annealing. 

Figure  6  shows  the  dislocation  substructure  in  a  Si -doped  GaAs  crystal 
after  it  was  deformed  15  pet  at  500°C  and  annealed  in  an  evacuated  and  sealed 
silica  tube  at  800°C  for  one  hour.  The  annealing  treatment  has  caused  a 
considerable  reduction  in  the  density  of  dislocations.  Dislocations  have  been  able 
to  react  during  the  annealing  treatment  to  produce  networks.  Another  feature 
worth  noting  is  the  occurance  of  many  small  prismatic  loops.  This  structure 
is  regarded  as  "recovered"  and  not  "recrystallized"  because  the  lack  of 
misorientation  between  regions  could  be  interpreted  as  grains. 

Figure  7  shows  the  dislocation  substructure  of  a  Si-doped  GaAs  crystal 
o 

deformed  20  pet  at  550  C  and  process  annealed  in  an  evacuated  and  sealed  silica 
tube  for  10  minutes  at  I100°C.  This  substructure  which  is  regarded  as 
recrystallization  consists  of  extensive  networks  of  dislocations  bounding  regions 
that  have  a  significant  misorientation  with  respect  to  one  another.  This  niis- 
orientation  is  evidenced  by  the  change  in  contrast  upon  crossing  a  dislocation 
network  which  can  be  seen  in  Fig.  7. 
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Fig.  6  Recovered  structure  of  Si-doped  GaAs  deformed 

15  pet  at  500°C  and  annealed  for  one  hour  at  800°C 


Fig.  7  Recrystallized  structure  of  Si-doped  GaAs,  deformed 

20  pet  at  550°C  and  annealed  for  ten  minutes  at  1  ]  00°C 


23 


All  thermomechanical  studies  to  date  have  been  carried  out  on  single 
crystals.  Recrystallization  of  a  initially  polycrystalline  material  by  cold 
working  and  process  annealing  would  have  been  somewhat  easier  because  in  the 
recrystatlization  process  grain  boundaries  are  known  to  be  the  preferred  sites 
for  nucleating  new  grains.  It  is  believed  that  recrystallized  GaAs  may  well 
show  improvements  ii  absorption  in  comparison  to  melt  grown  single  crystals. 
This  is  because  the  recrystallization  process  may  lead  to  a  redistribution  of 
impurities  along  grain  boundaries. 

Cold  working  and  process  annealing  experiments  are  continuing  and  will 
be  extended  to  the  Cr-doped  and  high  resistivity  undoped  materials.  The  effect 
of  recrystallization  on  the  absorption  coefficient  of  these  materials  will  be 
determined. 
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d.  1  Surface  and  Interface  IR  Absorption 


Clarence  C.  Crowell,  T.  Mangir 

During  this  quarter  we  have  continued  our  efforts  to  build  and  test  the 

MIS  device  which  is  described  in  detail  in  previous  reports.  This  device 

configuration  would  supply  us  with  information  about  surface  states  at  the  air- 

11  -2 

GaAs  interface,  (which  reportedto  have  density  of  the  order  of  10  cm 
for  fast  states  by  Pilkuhn*). 

To  date,  we-have  not  been  able  to  observe  the  variation  of  MIS  capacitance 

with  gate  voltage  and  frequency,  and  we  have  not  yet  obtained  desired  (~1000  A) 

2 

flatness  on  a  large  enough  sample  area  (~  1  cm  ).  Because  of  these  results, 
we  have  undertaken  a  closer,  more  detailed  study  of  the  effect  of  surface 
preparation. 

For  this  purpose  we  have  examined  the  surfaces  using  infrared  micro¬ 
scopy,  after  each  stage  of  preparation.  The  preparation  steps  are  as  follows: 

1.  GaAs  wafers  are  degreased  10  min.  in  warm  TCE  and  10  min.  in 
boiling  Methanol,  then  etched  in  a  1  H^O^:  3  H^SO^:  i  H^O  solution 
for  6-7  minutes.  Figure  1  shows  the  wafer  at  this  stage. 

This  step  removes  most  of  the  saw  damage. 

2.  The  wafers  are  then  mounted  on  a  polishing  block  and  polished  with 
Microlite  for  about  1  hr. 

3.  They  are  then  rinsed  with  neutralizing  soap  and  DI  water. 

4.  Finally,  the  samples  are  etched  in  1  H^O^:  3  H^SO^:  1  II^O  ^  rnin. 

Steps  2  to  4  are  repeated  to  achieve  a  better  flatness.  The  surface 
after  completing  step  3  for  the  first  time  appears  as  in  Fig.  2. 

Fig.  2;  The  definite  orientation  of  wires  on  the  surface  is  due  to  the 
saw  damaged  layer.  It  is  also  obvious  that  there  is  a  differently  absorbing 
region  which  appears  as  the  brighter  dot  on  the  picture. 

Fig.  3:  It  shows  the  surface  after  completing  step  4  for  the  first  time. 
Pattern  which  is  due  to  saw  damage  has  about  disappeared. 

After  repetition  of  steps  2-4  the  surfaces  appear  to  be  flat  to  the  order 
of  3000  under  multiple  beam  interferometer. 
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The  appeal ance  of  a  differently  absorbing  region  in  the  pictures  lead  us  to 
suspect  the  presence  of  some  inclusions,  or  precipitates  or  gross  crystal 
defects  in  the  material.  Such  nonuniformities  could  cause  the  flatness 
problems,  because  of  resulting  preferential  etxhing,  when  they  appear  on  the 
surface. 

Some  conducting  irregularly  shaped  objects  were  observed  on  the  GaAs 

surfaces  with  the  SEM  (cf  Fig.  5).  This  material  is  nominally  doped  to 
17  -3 

10  cm  and  was  chosen  for  the  MOS  studies.  The  exact  source  and  nature 
of  the  inclusions  is  not  yet  clear. 

At  present  we  are  trying  to  analyze  these  "impurities"  by  a  probe  tech¬ 
nique  and  get  some  new  material  to  work  with. 

REFERENCE: 

1.  Pilkuhn,  M.  H. ,  J.  Phys.  Chem.  Solids  2J5,  141  (1964). 


Dielectric  Constant  Measurement 

Clarence  R.  Crowell,  S.  Joshi 

In  the  course  of  this  research  program  it  became  necessary  to  know 
the  dielectric  constant  of  GaAs  accurately  to  predict  infrared  lbsses.  In  the 
literature  low  frequency  values  of  dielectric  constant  are  quoted  ranging 
from  10.  9  to  12.5.  The  frequency  and  temperature  dependence  of  dielectric 
constant  of  GaAs  would  also  be  of  interest  in  characterizing  the  material  as  a 
window  material  for  infrared  lasers.  For  these  reasons  an  experimental 
study  of  the  dielectric  constant  of  GaAs  was  undertaken. 

A  simple  way  to  measure  the  dielectric  constant  of  an  insulator  is  to 
measure  the  dielectric  capacitance  of  the  material  between  two  parallel  con¬ 
ducting  plates  which  act  as  electrodes.  In  a  practical  semiconductor  whose 
resistivity  is  finite,  the  measurement  frequency  should  be  high  enough  to  avoid 
effects  due  to  the  dielectric  relaxation  time  associated  with  free  carrier  re¬ 
distribution.  Also  in  order  to  avoid  the  effects  of  deep  impurities,  the  capa¬ 
citance  measuring  frequencies  should  be  in  excess  of  that  corresponding  to 
the  deep  level  response  time.  This  necessitates  that  the  material  be  of  as  high 
a  resistivity  and  as  little  compensated  as  possible.  Measurements  can  be  per- 
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formed  at  lower  temperature  at  which  carrier  concentration  and  deep  level 
response  time  will  be  reduced. 

Another  problem  associated  with  the  metal- semiconductor  system  is 
the  intrinsic  potential  barrier  associated  with  the  metal- semiconductor  contact. 
To  ensure  that  the  depletion  layers  under  the  contacts  extend  throughout  the 
sample,  the  Debye  screening  length  associated  with  the  carriers  must  be  much 
larger  than  the  sample  thickness,  i.  e.  the  effective  dielectric  capacitor  thick¬ 
ness  must  be  the  thickness  of  the  sample.  To  further  assist  in  ensuring  the 
above  condition  one  should  also  use  contacts  that  inject  neither  holes  nor 
electrons  into  the  GaAs.  Typical  Schottky  barriers  on  GaAs  should  satisfy 
this  criterion.  A  good  final  check  on  such  data  should  however  be  that  the 
final  capacitances  should  be  independent  of  applied  dc  bias. 

8 

Bearing  the  above  criteria  in  mind  a  GaAs  sample  of  about  10  ohm- cm 
resistivity  with  intrinsic  carrier  concentration  was  selected.  The  sample  was 
lapped  to  a  thickness  of  about  5  mils.  Gold  electrodes  were  then  evaporated 
on  both  sides  of  the  sample.  Preliminary  measurements  have  been  performed 
at  room  temperature.  The  sample  preparation  involved  two  deposition  cycles 
and  one  of  the  contacts  appeared  ohmic.  Consequently  the  capacitances 
measured  were  slightly  sensitive  to  bias. 

Further  measurements  need  to  be  performed  on  additional  samples  with 
rectifying  contacts  on  both  sides  and  if  necessary  at  lower  temperature. 
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d.  2  Study  of  Defects  in  II- VI  Compounds 

F.  A.  Kroger,  S.S.  Chern,  H.  R.  Vydyanath 


1. 

2. 

3. 

4. 


Experiments  on  the  self  diffusion  of  Cd  in  CdTe  •  In  at  high  indium 
contents  have  been  completed;  a  structure  in  the  isotherms  indicated 
by  earlier  experiments,  turns  out  to  be  real. 

The  results  are  being  analyzed  in  terms  of  a  detailed  defect  model. 


Species  to  be  included  are: 

V*'  ,  Cd!*,  V"  ,  Te.',  Te*.  In'  In*'  e'  .  h* 

Te  l  Cd  l  l  Cd  Cd 

The  introduction  of  doubly  ionized  In^  is  necessary  to  explain  an 
overshoot  effect  in  the  Hall  effect  first  tentatively  attributed  to  Inj '  . 
Though  In  of  a  higher  valency  is  in  general  not  found,  its  existence 


is  indicated  by  compounds  such  as  In  Te  . 

2  5 

Steps  are  being  taken  to  get  material  that  allows  us  to  study  electrical 


and  optical  properties  of  cooled  crystals. 

Further  study  of  the  function  of  Ag  in  CdS  has  shown  that  at  all  Cd 


pressures,  low  concentrations  of  Ag  act  as  donor;  at  low  p  higher 

U  Q 

ones  act  as  acceptors,  the  electron  concentration  concentration  as 
f([  Ag  ]  )  passing  through  a  maximum. 

At  still  higher  concentration,  the  Fermi  level,  stabilizes  approximately 
at  the  level  expected  in  undoped  CdS. 
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e.  1  Theoretical  Studies  of  Absorption  Mechanisms  in  IR  Window  Materials 

R.  W.  Hellwarth  and  M.  Mangir 

e.  l.a  Checks  of  Multi- phonon  Absorption  Theory 

_ 

During  this  quarter  we  carried  out  calculations  of  weighted  averages 
over  frequency  (or'moments1)  of  the  multi-phonon  absorption  in  a  rocksalt 
structure  crystal.  All  previous  calculations  of  multi-phonon  absorption  at  a 
given  frequency  have  required  drastic  approximations  of  unknown  error.  Our 
'moment1  calculations  however  are  exact  to  a  given  order  in  the  anharmonicity 
parameter,  and  have  printed  out  several  types  of  errors  in  existing  multi¬ 
phonon  absorption  theory.  For  example,  we  have  found  a)  that  absorption 

of  a  photon  leading  to  the  production  of  an  even  number  of  phonons  should  vanish 
in  the  limit  of  very  unequal  ionic  masses,  and  b)  the  actual  temperature  depend¬ 
ence  of  a  given  moment  is  different  than  that  which  would  have  been  predicted  for 
the  same  interaction  potential  by  previous  theories.  The  details  of  this  work 
and  itr  results  are  set  forth  in  the  accompanying  Appendix.  We  believe  the 
calcu. ation  of  moments  will  develop  into  a  valuable  tool  for  understanding  multi¬ 
phonon  absorption. 

e.  1.  b  Linear  and  Nonlinear  Electric  Moments 

This  quarter  we  continued  the  search  for  practical  experimental  methods 
for  determining  whether  or  not,  in  a  given  window  material,  multi-phonon  ab¬ 
sorption  is  mediated  mainly  or  in  part  by  a  nonlinear  relation  between  the  macro¬ 
scopic  eiec:iic  moment  per  unit  volume  (i.  e.  ,  polarization  density)  and  the  local 
ionic  displacements.  No  clearly  practical  methods  have  yet  been  found,  and 
the  search  continues.  In  this  regard  we  are  collaborating  with  Dr.  Crowell  on 
this  project  on  the  interpretation  and  re-determination  of  the  dc  dielectric 
constant  of  GaAs  and  other  window  materials.  There  is  great  disagreement  on 
this  quantity  in  the  literature  and  it  does  reflect  the  presence  or  absence  of  a 
nonlinear  moment. 

e.  1.  c  Multi-Phonon  Absorption  Calculations 

In  the  second  quarterly  report  of  this  project  we  gave  numerical  results 
for  the  multi-step  absorption  of  a  slightly  anharmonic  diatomic  molecule.  A 
lattice  of  these  non-interacting  molecules  was  argued  to  have  a  'multi-phonon' 
absorption  spectrum  which,  when  smoothed,  should  resemble  that  of  a  real 


lattice  of  interacting  units  cells.  In  fact  we  found  that  a  degree  of  molecular 
anharmonicity  which  was  consonant  with  that  known  to  exist  from  bulk  modulus 
and  other  date^  reproduced  well  the  exponential  fall-off  with  frequency  of  the 
wing  absorption  that  has  been  documented  experimentally  by  Deutsch.  At  the 
time  we  noticed  a  curiousity  of  the  anharmonic  molecule  absorption  which,  for 
the  record,  we  reproduce  below. 

When  the  exact  vibrational  line  spectrum  of  such  a  molecule  is  plotted, 
unsmoothed,  the  locus  of  the  line  peaks  gives  a  startingiy  realistic  plot  of  a 
typical  alkali  halide  absorption  spectrum,  from  an  asymmetrically  broadened 
fundamental  absorption  band  to  the  subsidiary  multi-phonon  peaks  that  are  in¬ 
creasingly  broadened  and  washed-out  at  higher  phonon  numbers.  A  typical  such 
molecular  absorption  spectrum  is  shown  in  figure  1  with  the  molecular  parameters 
adjusted  to  fit  the  measured  wing  abosrption  of  KC1  crystal  at  300K.  That  the 
resemblance  of  this  to  a  crystal  spectrum  should  exist  when  the  ionic  masses  are 
very  different  is  printed  out  in  Appendix  I.  (In  this  case  each  light  ion  essentially 
vibrates  independently  in  a  rigid  cage  of  heavy  ions).  However,  we  have  no  proof 
that  the  resemblence  is  more  than  accidental  for  a  crystal  like  KC1  whose  masses 
are  close. 


Caption  for  Figure  1 


Fig.  1.  A  plot  of  the  line  spectrum  of  a  diatomic  nolecule  inter¬ 
acting  via  a  Morse  potential  with  anharmonicity  parameter  V,  resonant 
frequency  U,TO'  and  dipole  matrix  element  adjusted  so  that  (3  from  U>= 
to  5  UtjQ  is  nearly  the  measured  absorption  coefficient  of  a  KCl  crystal  at 
300  K.  The  modified  absorption  coefficient  p  is  defined  as  4nNa"(U))/nc 
where  a"  is  the  imaginary  part  of  the  molecular  polarizability  at  U),  N  is  the 
number  density  of  KCl,  and  N  is  the  refractive  index  of  KCl. 
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Checks  of  Multi-Phonon  Absorption  Theory 


* 


Robert  Hellwarth 

Departments  of  Electrical  Engineering  and  Physics 
University  of  Southern  California 
University  Park 

Los  Angeles,  California  90007 


There  is  considerable  uncertainty  as  to  the  accuracy  of  calculations  to  date  of 
multi-phonon  absorption  in  infrared  window  materials.  Here  we  show  that  s-uch 
uncertainties-can  be  checked  by  calculating  various  weighted  integrals  over 
frequency  ("frequency  moments")  of  the  imaginary  part  x  "ito)  of  the  electrical 
susceptibility,  to  which  the  absorption  coefficient  £((».*)  cm  1  is  related.  We 
obtain  expressions  for  six  moments,  evaluating  a  moment  of  the  2-phonon 
contribution  explicitly  for  an  anharmonic  pair -interaction  model.  This  exact 
result  reveals  inaccuracies  in  the  dependence  of  absorption  on  ion  masses 
and  temperature  predicted  by  approximate  theories. 

Key  Words:  Crystal  absorption  theory,  frequency  moments,  infrared 
absorption,  infrared  susceptibility,  infrared  window  materials,  multi¬ 
phonon  absorption. 


I.  Introduction 

There  has  been  a  number  of  theoretical  efforts  to  estimate  the  intrinsic  infrared  (IR)  absorption 
coefficient  BfoJjcm*1  of  IR  window  materials  in  the  high-frequency  "tail"  of  the  fundamental 
lattice  absorption.  [1-7]1  This  absorption  is  presumably  "multi-phonon";  that  is.  an  absorbed 
photon  is  converted  to  m-phonons  so  as  to  conserve  energy  and  crystal  momentum,  and  obey 
certain  selection  rules.  [8, 9]  Numerical  predictions  from  the  m-phonon  theories  for  8  have  had 
to  rely  on  rather  drastic  approximations  that  leave  uncertain  the  magnitude  of  the  theoretical 
errors.  [1 -7]  Here  we  point  out  that  various  weighted  integrals  over  frequency  of  the  imaginary 
part  x"(«)  the  lattice  electric  susceptibility  can  be  evaluated  exactly,  or  at  least  to  known 
accuracy.  Since  x"  is  related  to  g,  the  predicted  values  of  these  "moments"  can  be  used  to 
check  the  accuracy  of  the  theories  of  g.  We  perform  such  a  check  of  2-phonon  absorption  in 
a  diatomic  cubic  crystal,  and  point  out  errors  in  existing  theoretical  results  for  the  mass  and 
temperature  dependence  of  this  absorption.  For  example,  we  find  that  a)  this  absorption 
should  be  small  when  the  ionic  masses  are  very  different,  and  b)  when  the  masses  are  equal, 
the  variation  of  the  2 -phonon  absorption,  in  going  from  low  to  high  temperatures,  is  really 
about  twice  that  predicted  by  previous  approximate  theories. 

In  Section  2,  we  review  the  quantum  expressions  for  B  and  x  (the  electric  susceptibility). 

In  Section  3,  we  define  the  moments  of  X(b')  and  obtain  the  quantum  formulae  for  six  of  them, 
in  terms  of  averages  over  various  derivatives  of  the  interionic  potential,  assuming  a  linear 
relation  to  exist  between  the  dipole  moment  operator  and  the  ionic  displacements  (i.e.,  a  "linear 
dipole  moment").  We  derive  as  an  example  in  Section  4  the  fifth  frequency  moment  of  the 
2-phonon  absorption  for  an  anharmonic  pair  potential.  Then  in  5  ,  we  compare  our  predicted 
mass  and  temperature  dependence  with  that  from  the  approximate  theories  of  g. 


* 

This  research  was  sponsored  by  Defense  Advanced  Research  Projects  Agency  Contract  No. 
F19628-72 -C-0275  monitored  by  Air  Force  Cambridge  Research  Labs. 


Figures  in  brackets  indicate  the  literature  references  at  the  end  of  this  paper. 
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2.  Quantum  Theory  of  Absorption  and  Susceptibility 

At  optical  frequencies  well  below  the  bandgap,  we  may  picture  a  window  material  as 
composed  of  ions  (labelled  by  greek  indices  a, having  masses  n\-j,  and  interacting  via  a 
potential  V  that  is  a  function  of  their  position  vectors  r  .  The  electronic  interactions  are 
absorbed  in  the  ion  dynamics  via  the  Born -Oppcnheimer  approximation.  This  .gives  an  effective 
Hamiltonian  H  for  a  subvolume  of  the  material  (where  L  is  much  less  than  an  optical  wave¬ 
length  but  much  greater  than  the  interionic  spacing)  of  the  form 

H  =  <» 

where  the  pi|(j  are  the  components  (i=x,  y,  z)  of  the  momentum  operator  for  the  ion  a.  The 
potential  V  does  not  include  long  range  interactions  between  ions  in  LJ  and  in  neighboring  -re¬ 
gions.  These  interactions  are  assumed  to  be  mediated  by  the  long  wavelength  components  of 
the  electric  field  E,  that  is,  by  the  macroscopic  electric  field,  the  transverse  part  of  which  is 
the  propagating  field  whose  absorption  coefficient  we  seek.  The  field  interacts  with  the  sub¬ 
volume  Lr  via  a  perturbation  -ME-|x  whore  \  is  the  electronic  susceptibility  with  ions 
fixed,  and  M,  the  component  of  the  electric  dipole  moment  in  along  E,  is  taken  to  be 

M  =  T  e  x  (2) 

a  o  a  '  ' 

where  ea  is  effective  charge  of  the  ion  a,  and  xa  is  the  component  of  the  displacement  from 
equilibrium  of  ion  o  along  the  direction  of  E.  We  neglect  here  the  terms  in  M  that  arc  not 
linear  in  the  x^  as  this  is  thought  to  be  accurate  for  most  polar  insulating  materials  and  has 
been  assumed  in  most  theories  of  multi-phonon  absorption. 

The  complex^electric  susceptibility  x-  ai  frequency  W  of  an  optically  isotropic  crystal  is 
defined  as  1 

X:(«)  -  L'3<M(W)>/E(U>)  (3) 

1  3 

Where  <M(W)>  is  the  expected  amplitude  of  the  electric  moment  in  L  oscillating  in  response  to 
the  macroscopic  field  of  amplitude  E(W).  The  absorption  coefficient  is  related  to  X  =  XC  +Xj  by 

g(u>)  Im  [1  +  4n  x(te)]£  (4a) 


4nX!l(q>)q; 


Here  x"  is  the  imaginary  part  of  Xj,  .c  is  the  velocity  of  light  in  a  vacuum  and  n  the  refractive 
index.  We  now  write  the  well-known  quantum  mechanical  formula  for  x"  *n  terms  of  the  ion-ion 
response  functions  f°r  which  we  will  derive  simple  moment  formulae: 


where 


X"(«)  =nL'  ra.e  SnB(*)caee 

Sa  e(W)  S  Zrn,  n(Pm*  P,i)<m  KJn><n'x  Jm> 


Here  |m>  is  an  eigenstate  of  the  Hamiltonian  II  for  the  subvolume  L3  having  energy  En1  and 
probability  l\n  of  occupation.  En-Em.  It  is  from  this  exact  formulation  (5)  and  (6)  for 

X"  that  previous  theories  of  6(W)  have  started.  We  proceed  instead  to  derive  from  (6) 
formulae  for  moments  of  S^p  and  thereby  for  x"- 

3.  Moments  of  lon-lon  Response  Functions1 
The  l  moment  of  the  ion-ion  response  function  S^^te).  as  defined  by 


rm  J  .  „ 


coth  (2kT')Sag  (w>*Xcven- 


is  seen  from  (6)  to  equal 

^  ^ ’ 1  7  m,  nPm<  En '  Ern  > ^  <  ,ri 'xa  t ur- 1 ' x  B 'm>  ^  ^ 

for  thermal  equilibrium  when  Pn/Pm~cxp  (E  -E  )/kT.  From  (8)  it  follows  directly  that 
S  „  ^  -  2<x  x  _>/ fl 


5  „  /m  (Thomas -Kuhn-Reich  sum  rule) 
oe  a 


I 


l 


Hd/varth  f>-  3  c~  & 


(2)  _ 
aB 

<xo  (3V/axp)>/(meM 

(H) 

(3)  . 
aB 

2 

<3  V/bx^ 3xn> /mm 
a  b  a  p 

(12) 

(4)  _ 
aB 

2<  (3  V/3xa)(3  V/bXp  )> 

(13) 

(5)  . 
aB 

Ea  .1 •= ,s2  v'*Vrial,a2v/sV  r  ial>/"’amema 

(14) 

by  repeated  use  of  the  identity 

(Em-En)<m|Aln>=<m|[H,A]jn>  (15) 


where  A  is  any  operator  and  H  is  the  Hamiltonian  operator  of  (1).  In  (14),  r.  is  the  displace¬ 
ment  of  ion  a  in  a  direction  i=x,  y,  z.  In  (9)  -  (14)  the  x,  are  displacements  parallel  to  the 
applied  electric  field. 


4.  Application  to  Pair  Interactions  With  Hard  Cores 


The  interaction  potential  V  of  ionic  crystals  is  known  to  be  fairly  well  represented  by  a 
superposition  of  pair  potentials 

v"*Ea.Sv,V  VVV  1161 

where  x  is  the  displacement  from  equilibrium  position  it  of  ion  a.  In  multi-phonon  absorption 
calculations,  the  potential  V  has  been  assumed  to  be  mad?  up  ofa  harmonic  part  V^,  quadratic 
in  the  ion  coordinates,  and  an  anharmonic  part  Van,  which  arises  from  the  hard-core  repulsion 
between  nearest  neighbors,  and  which  has  been  assumed  to  be  a  power  series 


V  =  I  a  l  ,  c  r  "/n! 
an  a,  0  r  =  3  n  a  3 


(17) 


in  the  relative  ionic  displacement  r  along  the  equilibrium  separation  vector  it  -R„.  (a  and  3 
are  summed  over  all  nearest  neighbor  pairs.)  Since  the  repulsive  part  of  v(r)  is  tho&ght  to  be 
roughly  of  the-  form  exp-r/p,  the  coefficients  c  are  thought  to  be  all  of  the  same  order  of 
magnitude,  and  of  the  order  of  a  dimensionlessanharmonicity  parameter  X  [10] : 


cn~X 


(18) 


where 


(19) 


M  is  the  reduced  mass  of  an  ion  pair,  and  is  the  transverse  optical  (or  reststrahl)  frequency. 
The  values  of  p  for  alkali  halide  crystals  have  been  determined  from  bulk  modulus  data  and 
yield  values  of  \  around  0.3. 


Given  some  natural  assumption  about  the  exact  values  of  c  (e.  g.  ,  c  =  X  x  constant),  the 

present  theories  for  the  absorption  8  are  not  able  to  proceed  with  known  accuracy  even  to 

calculate  two-phonon  absorption;  ttic  k-space  integrals  are  too  complicated.  However,  we  can 
proceed  to  find  the  moments  ( 9 ) - ( 1 4 )  to  some  desired  order  in  the  anha rmonicity  parameter. 

We  illustrate  this  here  by  finding  the  part  of  the  fifth  moment  that  is  of  order  X^  and  repre¬ 
sents  multi-phonon  absorption. 

(5)  .  2 

Ajp  examination  of  S  in  (14)  shows  that  the  only  X  terms  are  those  proportional  to  c-,c^ 

or  c^  .  (Here  we  have  used  C2  to  symbolize  the  harmonic  force  coefficients  in  analogy  wit*li 

the  i.iuo,  ~>onic  coefficients  cn  for  n>2  in  (17).)  An  examination  of  the  structure  of  perturba¬ 
tion  theory  shows  that  to  order  X^,  the  C2C4  term  describes  a  shift  and  broadening  of  the  one- 
phonon  absorption  while  the  corresponding  c 3^  term  gives  rise  only  to  multi-phonon  absorption. 
In  fa  ct  it  describes  the  rate  of  all  combinations-  of  two-phonon  absorption  and  emission  among 
all  phonon  branches.  This  is  evident  because  1:3  is  the  coefficient  of  all  products  of  three 
phonon  creation  and  annihilation  operators  in  the  perturbing  potential,  and  therefore  in  t ho  C3 
term  in  y  ",  These  terms  contribute  a  factor  to  the  absorption  matrix  element  that  describes 
the  annihilation  of  t ho  k-0,  TO  phonon  (created  by  a  photon)  along  with  the  annihilation  or 
creation  (in  any  combination)  of  two  other  phonons  from  any  branches.  We  now  calculate  this 
"2 -phonon"  contribution  ^ph  to  the  fifth  moment  S'^'  of 

Inserting  the  n  =  3  term  of  (17)  into  (14),  performing  the  indicated  diffe rentation,  and 
summing  over  ions  as  indicated  in  (5)  gives  for  a  rocksalt  structure  crystal 

■2 


M 


//.‘Av.v rill  P-^ ft (p 


S2ph<5,*Jo  dWto>SXM(W)2.ph0n=T'c3SV;,L‘:,Sli<(Xa*Xa-l,S 

-<(xa-l*xa){xa+l-xa)>M/ma  1* 


2  2-3.  -3. 


(20) 


Here,  x  j  ate  the  relative  x-displacements  of  the  nearest  neighbors  to  ion  a  in  opposite 
directions  along  the  crystal  x-axis;  e  is  the  magnitude  of  the  effective  charges*  ea  and  eg,  and 
M  is  the  reduced  ion  pair  mass.  The  thermal  averages  indicated  by  <  >  are  needed  here  only 
in  the  harmonic  approximation  (X  =0).  For  a  diatomic  crystal,  harmonic -crystal  quantum  theory 
gives  [10] 


<  >  .  -  »U-1(»..;Q)U(«,  8;Q)  coth 

a  9  harm  Q  2w  2hT 

Q  <1  P 


(21) 


where  the  Q  labels  all  the  normal  modes  of  the  crystal  in  volume  L' ,  and  Wq  are  their  i 
frequencies.  The  unitary  U-matrices  are  those  which  diagonalize  the  matrix  v 
deriving  from  the  expression  xy,  ap  p 


for  the  harmonic  part  of  the  potential. 


^h  ”  ^xy,a(ita/pvxy,  ap 


(22) 


5.  Comparison  With  Absorption  Theory 

There  exists  no  explicit  calculation  of  all  2 -phonon  processes  for  a  real  lattice  with  an 
explicit  with  which  we  might  compare  the  results  (20)  and  (21 )  directly.  However,  we  can 
obtain  some  interesting  checks  of  absorption  theory^examining  two  extreme  rocksalt  structure 
models  in  which  there  are  either  1)  very  unequal  ion  masses,  or  2)  equal  masses  and  nearest 
neighbor  interactions. 


51  ma  +  l>>ma 

When  the  masses  of  the  ions  in  a  rocksalt  cubic  crystal  are  very  different,  then  the  light 
ions  (a)  move  in  a  nearly  rigid 
This  with  the  limit  ma-*p  cause 
2-phonon  absorption  vanishes  when  m^  j»m^.  This  is  reasonable,  because  the  light  ion  is 
moving  effectively  in  a  potential  having  inversion  symmetry,  and  its  states  have  even  or  odd 
parity.  Since  dipole  transitions  are  not  allowed  between  states  of  the  same  parity,  no  even 
number  phonon  absorption  processes  will  exist  in  this  limit.  No  hint  of  this  depression  of  the 
even-phonon  number  absorption  is  found  in  present  approximate  theories,  although  the  effect 
can  be  seen  in  the  relative  experimental  heights  of  the  2-phonon  resonance  peaks  in  different 
alkali  halides.  [11]  This  failure  probably  signals  factors -of -two  errors  in  the  high-frequency 
tail  absorption  of  crystals  like  LiF  whose  ionic  masses  are  quite  different. 


cage  of  the  surrounding  heavy  ions,  making  x^j-xssx^  i*xa’ 
the  two  averages  in  (18)  to  cancel  exactly,  indicatinguthat 


5.2  ma 


ma+j-2(iand  Nearest  Neighbor  Forces  Only 


When  both  ion  masses  arc  equal  and  the  harmonic  forces  are  between  nearest  neighbors  in 
a  cubic  diatomic  crystal,  the  normal  lattice  modes  consist  of  three  degenerate  branches  whose 
ion  motions  are  polarized  along  x,  y  and  z.  Since  we  need  only  correlations  between  x-dis¬ 
placements  in  (20),  the  sum  in  (21)  is  over  one  branch  only.  The  normal  mode  coordinates  for 
this  brand,,  as  embodied  in  the  unitary  U-matrices,  arc  the  familiar  plane  waves 


U(x,  a ;Q)  =  (^)*  expi-2-^  (Qxax+Qyay+Qzaz) 

where  the  modes  are  labelled  by  the  integers  Q^,  Q  ,  :  0,i  1 ,  .  .  i  and  the  ions  by  the 

integer  coordinates  ax.tly,az  ~  1,2,...  I, /a.  The  IcYigtlY  of  tie  cube  cllgo  of  the  crystal  sub¬ 
volume  is  I.  and  the  interionic  spacing  is  a.  Diagonalizing  the  harmonic  potential  with  (23) 
gives  the  eigenvalues  of  the  potential  matrix  and  the  dispersion  relation  for  the  x-branch: 


(23) 


W  =  W  sin  (nQ  a/L), 
U  o  x 


(24) 


where  i 

r:  (2 C/m** 


(25) 


is  the  "TO"  or  fundamental  (k-0)  optical  phonon  absorption  frequency  for  the  rocksalt  structure 
lattice  with  equal  masses  (m-2|j)  and  nearest-neighbor  force  constant  C. 


He/lwsi'n  p-5“ ot  fc> 


Substituting  (23)  and  (24)  in  (21)  gives  for  the  2-phonon  absorption  moment  of  (20) 


,2  2. 

o  (5).  Zc3  e 
»-»_L  -  — i—j -  I(Y) 


2ph  4  . 

r  li  a  to 


where 

The  temperature  parameter  y  is 

At  high  temperatures 
and  at  low  temperatures 


o 


I(V)  s  JQ  du  u(l-u  coth  yu. 
.  huo 

Y=  2kT  ‘ 


I(yH  ^7*  y -*0 


I(YH  3  .  Y  • 


(26) 

(27' 

(28) 

(29) 

(30) 


For  fixed  B  rc^.e,  and  a,  the  independent-anharmonic -molecule  models  [1-4]  predicted  for  this, 
or  any  other  type  of  lattice,  a  2-phonon  temperature  factor  2n  +  1,  where  n  is  the  Planck 
function  [e^Y-l]**.  Therefore  these  very  approximate  theories  predict 


s  <5>(T.  0, 

~{$) -  =  y 

S2ph  <T~> 


(31) 


rather  than  the  exact  result  for  this  lattice  model  of  4y/3n  from  (29)  and  (30).  A  plot  of  I(y) 
between  its  asymptotes  is  shown  in  figure  1  together  with  the  2n+l  function  normalized  tc 
approach  I(y)  at  high  temperatures. 


It  is  well-known  that  this  model  having  only  nearest  neighbor  forces  does  not  accurately 
describe  any  real  rocksalt-structuro  ionic  crystals.  However,  we  feel  that  checking  a  theory 
by  applying  it  to  this  model  does  give  a  reasonable  estimate  of  its  accuracy  for  a  more 
realistic  model  of  a  crystal.  That  is,  the  fact  that  the  ratio  (31)  from  approximate  theories  of 
absorption  is  2.4  times  the  exact  result  for  this  (albeit  poor)  model  is  an  indication  of  the  kinds 
of  errors  for  the  moments  of  6  that  the  independent -molecule  models  may  yield  in  the  2-phonon 
region.  It  is  also  likely  that  the  errors  are  even  larger  at  higher  frequencies  (higher  phonon 
numbers). 

In  summary,  we  have  shown  how  calculations  of  frequency  moments  of  t ho  imaginary  part 
X"(w)  of  the  electric  susceptibility  can  be  performed  with  much  greater  amuaci  than  can 
calculations  of  \"(u:)  itseif,  and  how  the  moment  calculations  can  be  used  to  cheek  the  more 
detailed  calculations  of  multi-phonon  absorption  fi(u;)  as  a  function  of  frequency.  Many  other 
checks  and  relations  may  be  derived  from  the  moment  expressions  (9)  to  (14)  than  we  have 
mentioned  here. 
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7.  Figures 


Figure  1.  Plot  of  (he  temperature  factor  I(y)  calculated  for  the  fifth  frequency  moment 
of  all  2-phonon  contributions  to  (lie  imaginary  part  v'  '  of  the  susceptibility 
of  a  roc k sal t- st rue lu re  crystal  with  nearest  neighbor  interactions.  The 
factor  2n  1 1  predicted  from  the  independent-molecule  approximation  foi  the 
same  moment  is  shown  for  comparison,  normalized  arbitrarily  to  approach 
the  exact  result  at  high  tempt  rate  res .  (\  ’  fiU^/ZkT) 
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f.  1  Wavelength  Dependent  Calorimetry 

W.  H.  Steier,  S.  T.K.  Nieh,  R.  Joiner 

NEW  SAMPLES 

In  accordance  with  our  plan  of  locating  the  highest  purity  melt  grown 
GaAs  samples  grown  in  the  past  5  to  10  years,  we  have  requested  and  received 
samples  from  IBM  and  Bell  Telephone  Labs.  The  IBM  samples  were  prepared 
by  Jerry  Woodall  using  the  horizontal  bridgeman  technique.  Before  growth 
oxygen  was  purposely  introduced  into  the  boat  and  after  growth  the  sample  was 
remelted.  (Hence  the  name  WORM-1 ,  With  Oxygen  ReMelt. )  Particular  care 
was  given  to  purity  of  source  material  and  avoidance  of  contamination  during 
growth.  The  result  of  our  calorimetric  measurement  of  WORM-1  is  almost 
identical  to  that  of  WA1000  reported  in  our  last  quarterly  report. 

The  Bell  Telephone  Lab's  samples  were  prepared  using  the  float 
zone  technique  developed  by  Professor  James  Whelan.  Table  1  contains  in¬ 
formation  of  different  samples.  Preliminary  calorimetric  measurements  show 


that  the  residual  optical  loss  of  these  samples  are  similar  to  that  of  WA  1000. 
More  detailed  measurements  are  planned  to  confirm  this. 

Sample 

No. 

2 

cm 

V  sec 

Room  Temp.  Mobility 

ft -cm 
Resistivity 

Comment 

44- 3A 

5550 

4.  6x10  8 

carbon  doped 

54-4A 

5700 

2.  5xl08 

slow  pass 

55-1A 

6000 

1.  4x10 8 

slow  pass 

46-5A 

Inhomogeneous 

Inhomogeneous 

high  As  pressure 

Table  1  -  Bell  Telephone  Lab  Float  Zone  Samples 
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Conclusions 


After  having  measured  the  highest  purity  and  best  crystalline  quality  melt  grown 
GaAs  crystals  available  anywhere,  we  are  forced  to  conclude  that  there  is  a 
definite  lower  limit  on  the  absorption  coefficient  of  about  0.008  cm*'  at  10.  6n. 

The  exact  nature  of  the  basic  mechanism  responsible  for  this  loss  is  not  clear 
at  this  point.  However,  we  hope  tc  have  a  better  idea  after  we  have  obtained  the 
wavelength  dependent  absorption  data  as  a  function  of  temperature.  From  the 
fact  that  high  pGaAs  needs  heavy  compensation  we  can  speculate  that  impurity 
plays  an  important  role.  Localized  nonstoichiometry  could  be  another  possible 
mechanism. 

Dr.  Don  Stierwalt  of  NELC,  San  Diego  recently  kindly  consented  to  make  emis- 
sitivity  measurements  on  a  Bell  and  Howell  sample  EMC  6012B  which  has  al¬ 
ready  been  measured  by  us  calorimetrically.  Figure  1  shows  his  results  with 
our  data  superimposed  upon  it.  The  emissitivity  measurements  were  made  at 
373°K  and  the  calorimetry  measurements  at  300°K.  Other  measurements  near 
273°K  at  NELC  are  very  similar.  The  agreement  between  the  data  is  very  good 
and  provides  an  independent  check  on  both  measuring  techniques.  The  structure 
is  identical  except  for  a  shift  in  wavelength.  This  could  be  a  temperature  effect, 
but  is  probably  due  to  a  wavelength  calibration  error.  The  USC  wavelength 
calibration  is  believed  to  be  correct  since  it  agrees  with  the  known  output  spectra 
of  the  CO^  laser.  This  agreement  is  a  very  encouraging  result  for  both  the 
NELC  and  the  USC  groups. 

Computer  Simulation  of  Heat  Flow  in  Calorimetry 

The  two-dimensional  heat  program  TEMPS  is  now  working.  An  initial  run  was 
made  using  insulated  boundary  conditions  and  actual  laboratory  sample  para¬ 
meters.  The  computed  time  rate  of  change  of  the  surface  temperature  agreed 
within  5%  with  the  results  obtained  from  the  theory  used  in  the  calorimeter 
measurements. 

The  three-dimensional  heat  program  HEATING. 3  has  been  received  and  is 
being  made  ready  for  use. 
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Conclusions 


After  having  measured  the  highest  purity  and  best  crystalline  quality  melt  grown 
GaAs  crystals  available  anywhere,  we  are  forced  to  conclude  that  there  is  a 
definite  lower  limit  on  the  absorption  coefficient  of  about  0.008  cm  '  at  10.  6|i. 

The  exact  nature  of  the  basic  mechanism  responsible  for  this  loss  is  not  clear 
at  this  point.  However,  we  hope  to  nave  a  better  idee,  after  we  have  obtained  the 
wavelength  dependent  absorption  data  as  a  function  cf  temperature.  From  the 
fact  that  high  pGaAs  needs  heavy  compensation  we  can  speculate  that  impurity 
plays  an  important  role.  Localized  nonstoichiometry  could  be  another  possible 
mechanism. 

Dr.  Don  Stierwalt  of  NELC,  San  Diego  recently  kindly  consented  to  make  emis- 
sitivity  measurements  on  a  Bell  and  Howell  sample  EMC  6012B  which  has  al¬ 
ready  been  measured  by  us  calorimetrically.  Figure  l  shows  his  results  with 
our  data  superimposed  upon  it.  The  emissitivity  measurements  were  made  at 
373°K  and  the  calorimetry  measurements  at  300°K.  Other  measurements  near 
273°K  at  NELC  are  very  similar.  The  agreement  between  the  data  is  very  good 
and  provides  an  independent  check  on  both  measuring  techniques.  The  structure 
is  identical  except  for  a  shift  in  wavelength.  This  could  be  a  temperature  effect, 
but  is  probably  due  to  a  wavelength  calibration  error.  The  USC  wavelength 
calibration  is  believed  to  be  correct  since  it  agrees  with  the  known  output  sp'ectra 
of  the  CO2  laser.  This  agreement  is  a  very  encouraging  result  for  both  the 
NELC  and  the  USC  groups. 

Computer  Simulation  of  Heat  Flow  in  Calorimetry 

The  two-dimensional  heat  program  TEMPS  is  now  working.  An  initial  run  was 
made  using  insulated  boundary  conditions  and  actual  laboratory  sample  para¬ 
meters.  The  computed  time  rate  of  change  of  the  surface  temperature  agreed 
within  5%  with  the  results  obtained  from  the  theory  used  in  the  calorimeter 
measurements. 

The  three-dimensional  heat  program  HEATING  3  has  been  received  and  is 
being  made  ready  for  use. 
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ABSORPTION  COEFFICIENT 


f.  2  Alkali  Halide  Surface  Studies  with  Acoustic  Probe  Techniques 
J.  H.  Parks,  T.  Colbert 

A  previous  analysis ^  indicated  the  sensitivity  and  utility  of 
acoustic  surface  wave  propagation  for  the  study  of  alkali  halide  surface 
properties.  The  primary  emphasis  during  this  quarter  has  been  to 
develop  these  acoustic  probe  techniques.  This  effort  has  concentrated 
on  achieving  surface  wave  propagation  on  KC1  surfaces,  and  beginning 
experiments  to  study  the  details  of  acoustic  detection  of  infrared  absorption 
at  Y-cut  Quartz  surfaces. 

I  Acoustic  Wave  Propagation  on  Alkali  Halide  Surfaces 

The  interdigital  transducer  used  to  excite  surface  acoustic  waves 
relies  on  a  piezoelectric  coupling  at  the  material  surface  on  which  it  is 
placed.  For  this  reason,  it  is  not  possible  to  generate  surface  waves  on 
the  centrosymmetric  alkali  halide  surface  by  direct  application  of  these 
transducers.  Several  methods  can  be  applied  to  achieve  wave  propagation 
an  non  piezoelectric  surface  including 

(a)  bulk  wave  to  surface  wave  conversion, 

(b)  growth  of  a  piezoelectric  crystal  film  (eg.  sodium  chlorate,  lithium 
iodate,  aluminum  nitrite)  on  the  material  surface  which  then  serves 
as  a  base  for  the  transducer, 

(c)  sputtering  a  thin  film  of  zinc  oxide  on  the  non  piezoelectric  surface 
which  becomes  oriented  during  deposition  and  then  serves  as  a 
transducer  base, 

(d)  coupling  surface  acoustic  waves  from  a  piezoelectric  surface  such 
as  quartz  through  a  fluid  film  interface  to  the  surface  of  the  non 
piezoelectric  material. 

Our  initial  attempt  will  be  to  excite  surface  propagation  on  alkali 
halide  surfaces  via  fluid  coupling  from  Y-cut  Quartz.  This  method  has 
been  chosen  because  it  demands  the  least  sample  preparation  and  handling. 

This  is  desirable  since  we  plan  to  apply  these  techniques  to  study  ultra 
high  purity  alkali  halide  surfaces. 

A  device  to  study  the  propagation  behavior  of  acoustic  waves  which 
are  fluid  coupled  to  alkali  halide  samples  is  shown  in  Figure  la.  The 
precision  micrometer  head  allows  two  acoustic  transducers  on  individual 
Y-cut  Quartz  crystals  to  be  brought  into  contact  with  a  fluid  interface 
between  the  Quartz  and  sample  surfaces.  In  Figure  lb  the  mounting  of  the 
Quartz  and  sample  is  indicated  to  clarify  the  geometry  of  the  fluid  coupling 
interface.  The  Quid  layer  is  held  between  the  Quartz  and  sample  by  capillary 
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and  does  not  extend  further  on  the  sample  surface.  As  shown  in  Figure  lc, 
the  acoustic  wave  initially  propagates  within  this  fluid  layer  and  on  the  surfaces  1 
of  the  Quartz  and  a  ample.  This  wave  will  continue  propagating  beyond  the 
fluid  layer  as  a  surface  acoustic  wave  on  the  sample.  This  fluid  coupling 

( 

technique  has  been  studied  experimentally  in  References  (2)  and  (3)  and  the  i 

propagation  of  layered  waves  is  discussed  in  References  (4)  and  (5)  and  more 
recently  in  Reference  (6). 

We  are  beginning  an  experimental  study  of  surface  wave  propagation 
on  KC1  and  NaCl  samples.  Initially  the  propagation  behavior  will  be 
investigated  using  purchased  material  which  has  been  polished  to  reduce 
surface  irregularities  to  less  than  ~  l|i  size.  Polishing  techniques  are 
being  developed  similar  to  those  at  NRL  to  allow  us  to  prepare  samples 
grown  here  at  USC.  These  experiments  will  provide  the  following  information: 

(a)  The  effects  of  disspersion  and  acoustic  mode  conversion  introduced 
by  layered  structure,  including  the  wave  velocity  as  a  function  of 
fluid  layer  thickness, 

(b)  acoustic  frequency  limitations  caused  by  dispersion  and  also  by 
attenuation  of  the  material  surface, 

(c)  the  stability  of  the  fluid  layer  thickness. 

The  last  point  (c)  is  essential  for  the  measurement  of  small, 
acoustic  wave  phase  variations.  ^  Mechanical  vibrations  and  evaporation 
of  the  fluid  which  change  the  fluid  thickness  will  introduce  variations  in 
the  phase.  For  example,  a  phase  variation  of  ~10  m  rad  corresponds  to  a 
change  in  the  fluid  layer  thickness  the  order  of  0.  5u.  However,  the 
phase  change  signal  induced  by  sample  absorption  will  have  a  time  dependence 
which  is  significantly  different  from  that  of  the  layer  thickness  variation 
and  will  allow  these  effects  to  be  separated.  Another  possibility  for 
eliminating  this  uncertainty  involves  the  use  of  two  parallel  acoustic  waves 
propagating  through  the  fluid  layer  structure.  Any  fluid  thickness  variation 
will  affect  both  waves  similarly  and  allow  a  null  detection  technique  to 
eliminate  these  effects.  The  separation  of  these  waves  by  a  distance  of  ~  1cm 
will  isolate  these  waves  from  the  coupling  effects  of  thermal  diffusion  in  the 
KC1,  NaCl  surface  absorption  experiments. 

II  Acoustic  Detection  of  10,  6u  Radiative  Absorption 

Initial  experiments  have  been  performed  to  detect  infrared  absorption 
using  acoustic  wave  techniques  described  in  a  previous  report.^ 


These 


experiments  are  providing  information  on  the  possibility  of  using  these 
techniques  to  develop  a  new  type  of  infrared  detector,  and  also  to  develop 
a  method  of  measuring  surface  absorption  which  will  be  applied  to  alkali 
halide  materials.  In  these  experiments  10.  6p  CC>2  laser  radiation  was 
detected  via  surface  wave  propagation  on  Y-cut  Quartz  substrates.  The 
experimental  setup  consisted  of  a  CC^  laser  which  passed  through  a  camera 
shutter  to  generate  pulses  of  10  m  sec  to  100  m  sec  duration.  An  internal 
aperature  was  used  to  control  beam  power  and  improve  spatial  mode  quality. 
To  retain  beam  uniformity  this  radiation  was  not  focused,  and  the  detector 
shown  in  Figure  2  surface  was  limited  by  a  3  mm  diameter  aperature.  The 
interdigital  transducers  generated  a  120  MHZ  surface  wave  of  wavelength 
^=25  which  travelled  through  a  pathlength  of  2  cm  (of  which  only  3  mm 
was  exposed  to  infrared  radiation).  The  spatial  distribution  of  the  laser 
intensity  varied  with  power  and  limited  the  accuracy  of  determining  the 

actual  10.  6|i  intensity  probed  by  the  ~  1  mm  wide  acoustic  beam  to  ~±25%. 

2 

The  mismatch  of  laser  irradiated  area  (~  7  mm  )  and  acoustic  probed  area 
3 

(~3  mm  )  reduced  the  effective  infrared  intensity  depected  by  a  factor  of 
~  0.40. 

Figures  3a  and  3b  show  representative  data  traces  of  detected 
radiation  for  various  laser  powers  and  pulsewidths.  The  detected  laser 
power  P,  radiative  pulsewidth  t^,  the  maximum  phase  variation  are 

indicated  below  for  these  data  traces.  The  acoustic  wave  phase  variation 
was  measured  using  a  HP  8405A  Vector  Voltmeter  and  displayed  on  a 
Tek  549  Storage  Oscilloscope,  The  lKHz  bandwidth  of  the  Vector  Voltmeter 
limited  the  detector  response  and  accordingly  the  initial  data  includes  only 
10.  6  p  pulses  having  10  m  sec. 

Figure  3a  shows  a  small  phase  variation  of  A^j^O.7  mrad  induced 
by  a  10.  6p  radiation  with  a  timescale  of  50  m  3ec  per  division.  The 
multiple  exposure  shown  in  Figure  3b  represents  a  set  of  detected  pulses 
each  having  P=-0.  2  watt  and  different  pulsewidths  of  Tp=20,  28,  47  and  90 
m  sec  respectively.  Each  detector  signal  shows  a  linear  ramp  which  is 
characteristic  of  the  temperature  increase  following  a  square  radiation 
pulsewidth.  The  effects  of  radial  heat  diffusion  can  be  observed  clearly 
for  the  90  m  sec  pulse  as  a  slight  curviture  in  the  increasing  signal  trace. 
The  oscilloscope  trace  timescale  is  50  m  sec  per  division  and  the  phase 
variation  scale  is  1.5  mrad  per  division. 

The  detection  of  long  radiative  pulsewidths  used  in  these  initial 
experiments  cannot  be  analyzed  by  heat  equation  solutions  which  neglect 
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thermal  diffusion  terms.  Although  radial  diffusion  cut  of  the  acoustic  beam 

path  is  «  100  msec,  the  heat  diffusion  perpendicular  to  the  surface  requires 

only  ~4psec  for  heat  to  diffuse  out  of  the  depth,  d^ZSpt.  This  diffusion  into 

the  crystal  bulk  can  reduce  the  sensitivity  by  a  factor  of  ~  50  for  a  10  m  sec 

radiation  pulse;  however  short  pulses  <_  sec  will  be  detected  with  full 

sensitivity.  Analysis  of  the  infrared  detection  data  is  presently  underway 

using  numerical  solutions  of  the  heat  equation  generated  by  Heating  3: 

(7) 

An  IBM  Heat  Conduction  Program' 
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Figure  lb.  Fluid  coupling  arrangement  of  transducers 
and  sample 


KC1  Sample  Interdigital  Transducer 
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Figure  lc.  Fluid  coupling  schematic  showing  acoustic 
wave  propagation 


Figure  2.  The  acoustic  surface  wave  propagates  between  the 
■  two  upper  transducers  separated  by  2  cm.  These 
transducers  are  on  a  1.  5  inch  square  substrate  of 
Y-cut  crystal  quartz  which  is  bonded  to  an  electrode 
connector  holder. 


i 


Figure  3a.  Acoustic  wave  phase  signal  induced  by 

absorption  of  a  50  mwatt  pulse  having 
a  pulsewidth  of  20  msec. 


Figure  3b.  Multiple  exposure  showing  acoustic  wave 
phase  signals  inducer!  by  the  absorption 
of  .2  watt  pulses  having  pulsewidths  of 
20,  28,  47,  and  90  msec,  respectively. 
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g.  1  Characterization  of  Optical  Performance  of  IR  Window  Systems 


Martin  Flannery,  J.  H.  Marburger 

In  the  last  quarterly  report  we  described  an  approximate  method  for 
treating  window  birefringence  from  a  narrow  laser  beam  if  the  beam  could 
be  approximated  by  the  superposition  of  several  circularly  symmetric  beams 
clustered  near  the  window's  center.  The  method  is  simply  to  consider  each 
circular  piece  of  the  beam  to  be  centrally  located  in  a  disc  and  to  solve  for 
the  temperature  distribution  and  resulting  stresses  on  that  disc,  and  finally 
to  superpose  the  solutions  for  all  pieces  of  the  beam.  The  object  of  the 
approximation  was  to  find  relatively  accurate  solutions  of  slightly  asymmetric 
beams  without  the  necessity  of  using  a  general  thermoelastic  computer  p-ogram. 

It  was  decided  to  use  this  approximation  to  include  birefringence  in  a 
further  investigation  of  the  beam  slewing  reported  for  a  KC1  window  by 
Dr.  L.  Skolnik  of  AFCRL.  A  central  scan  of  the  laser  beam  used  in  these 
experiments  was  provided  by  Dr.  Skolnik  as  described  in  the  2nd  Quarterly 
Report. 


To  match  this  central  scan  we  developed  the  following  intensity  profile: 


(1) 

(2) 


I  (x.y)/Io=  i(r)  +  i(s) 


i(r)  =  A(e 


-OJr 


2 

■  iDa  )  4 


Br2(C-r)4  ;  r<C 
0  :  r>C 


(3)  i(s)  = 


-U)  s 

Cj  (e  V 


-U)  b 


(4)  s  =  (x-xo)2  4  y2 

where  r=0  is  the  window  center,  r=a.  The  window  edge  and  i(r)  the  central  beam. 

For  the  secondary  beam,  i{s),  centered  at  r=  xq,  we  will  assume  Cj,<<:A  and  that 

it  covers  a  disc  of  radius  b  =  a+xQ*  Both  pieces  of  the  beam  were  adjusted  to 

have  zero  intensity  at  the  edges  of  their  respective  discs  so  the  heat  equation  for 
each  disc  could  be  solved  with  the  boundary  conditions  T  (edges)  =  0. 

In  what  follows  we  will  generally  only  work  out  the  solutions  for  i(r)  since 

the  solutions  for  i(s)  can  be  gotten  from  them  by  a  suitable  change  of  parameters. 
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Also  all  distances  are  assumed  to  he  normalized  to  d,  which  is  usually  taken 
to  be  the  diameter  of  the  aperture  used  in  the  computer  diffraction  program 
(and  which  may  be  taken  as  d  =  a,  but  d  <a  is  often  chosen  to  save  computer 
time  on  very  narrow  beams). 

I  THE  HEAT  EQUATION 

To  get  sufficient  characterization  of  the  window  performance  it  seems 
sufficient  to  solve  the  heat  equation  for  short  times,  assuming  the  diffusion 
term  is  zero,  and  for  the  steady- state  case.  For  simplicity,  we  will  use  the 
boundary  conditions  T  (edges)  =0. 

2 

In  terms  of  the  normalized  time  T=  kt/d  ,  where  k  is  the  thermal  diffusion 
coefficient  and  with  the  temperature  normalized  to  Tq  =  3l  /k^»  the  general 
heat  euqation  is: 

(5)  dtT  =  ?2  T  ,  Q/kthTo 

For  short  times  we  have  the  solution 

161  T  T  = T 1  '*•*>  no 

th  o 

The  fact  that  this  solution  is  linear  in  time  will  result  in  the  whole  phase 
delay  through  the  window  being  linear  in  time  and  allow  us  to  calculate  many  time 
steps  with  relatively  little  work. 

For  the  steady  state  case  we  must  solve  each  circular  beam  separately  to 
preserve  its  circular  symmetry.  With  the  aid  of  the  Greens  function  for  a 
circular  symmetric  heat  distribution, 

(7)  G(r|r^)  =  -4tt  i.nr  ;  r>r^ 

=  -4TT  SL  nr^  ;  r<  r^ 

we  find  the  steady- state  solution  is: 

(8)  T(r)  =  —  G(r  |  r1)  i  (r1)  r^r1 

4tt  £ 

and  similarly  for  i(s).  There  is  no  difficulty  in  integrating  (8)  for  the  intensities 
(2)  and  (3),  but  the  polynomial  in  (2)  does  produce  a  lot  of  terms  so  we  won't  write 
the  solution  here.  The  Gaussian  terms  in  (2)  and  (3)  produce  a  term  which  is  most 
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easily  written  in  terms  of  the  exponential  integral  function. 

II  OPTICAL  ABERRATION 

For  a  window  of  thickness  h  and  refractive  index  n,  the  optical  path  length  is 
•t=  nh.  When  the  window  is  heated  or  strained  the  change  in  the  optical  path 
length  to  first  order  is: 

(9)  A  i>-  (n,Q-l)Ah  -!  h  An 

where  the  factor  (n  -1)  accounts  for  the  change  in  path  length  on  expansion  as  the 
new  path  length  through  the  window  minus  the  old  path  length  through  the  air. 

There  are  a  large  number  of  effects  that  can  produce  changes  in  h  and  nQ 
but  we  will  consider  only  the  following: 

(A)  Free  thermal  expansion  in  the  direction  perpendicular  to  the  plane  of  the 
window; 

(10)  Ah  =  ah  T  T 

a  o  t 

where  a  is  the  linear  expansion  coefficient  and  Tfc  is  the  total  temperature 

(B)  Bulging  in  the  direction  perpendicular  to  the  window  due  to  the  partially 
constrained  expansion  in  the  plane  of  the  window  is  given  by: 

(11)  ihb  =  -vhaTo(Orr<  oM) 

which  follows  from  Hook's  law,  and  where  the  stresses  have  been  normalized  to 
aE  Tq,  E  is  Young's  modulus,  v  Poisson's  ratio.  For  an  isotropic  cylindrically 
symmetric  temperature  distribution  these  normalized  stresses  are: 


(12) 
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Using  these  expressions  in  equation  (11)  and  ignoring  a  constant  term  we  get 
the  constrained  bulging  as: 


(16)  Ah.  =  vh  a  T  T. 

D  O  t 

if  we  again  include  all  the  beam  components. 

The  refractive  index  will  change  because  of  its  normal  temperature  depend¬ 
ence  and  the  stress-optic  effect.  Using  Hook's  Law  we  can  write  the  incremental 
indicatrix  (in  terms  of  the  stresses)  for  a  cylindrically  symmetric  window  as: 

(17)  ABrr  =  n;2-„;2  =  <pu+2p12)  aToT  '  b|(0rr4  b^j 

(18)  AB39=  ne  -no  =  (put2p12)  AToT  +  bjla8e+  bJ_CTrr 

where  the  p's  are  the  strain-optic  coefficients  and  the  b's  are  stress-optic 
coefficients. 


(C) 

(19) 


The  normal  thermal  changes  in  the  refractive  index  are  usualy  expressed 
in  terms  of  the  experimental  value 


dn 

dT 


o  =0 


3 

n 


'  2  ({\r2  pi2* 


a 


so  the  isotope  terms  in  equations  (17)  and  (18)  can  be  included  here: 


(20)  An 


o=0 


dn 

dT 


a=0 


T  T. 
o  t 


(D)  The  remaining  anisotropic  terms  in  equations  (17)  and  (18)  are  the  stress- 

birefringence  effect.  There  will  be  a  set  of  these  equations  for  each  circular 
piece  of  the  beam,  but  in  order  to  get  the  total  incremental  indicatrix  we 
must  transform  them  to  the  same  coordinate  system,  which  we  choose  to 
be  the  cartesian  xy- system. 

The  anisotropic  part  of  the  refractive  indices  for  the  two  principle  waves 
cannot  be  found  by  diagonalizing  the  total  indicatrix  at  each  point  in  the  window. 

If  the  eigenvalues  are  AB.,  then  using  the  relationship 
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(?l)  AB.  =  n."2  -  n  "2  ^  2(n.-n  )/n^ 

11  o  1  o  o 

we  can  write  the  anisotropic  increments  in  the  refractive  index  as: 

(22)  An  ,  a  .  n  3  AB.  ,/2 

1, 2  o  1, 2 

The  total  path  length  changes  can  now  be  written  as 

(23)  ^2  =  (no-D  <Aha+  Ahfa)  4  hAn0=()  +  hAn^2 


from  equations  (9),  (10),  (16),  (20)  and  (22).  Theisotropic  part  is  characterized 
by  the  single  constant  G: 


(24)  Li 


[ 


hT  \  (n  -1)  (1+v)  a  + 


dn  1 

dT  0-0  J 


T  T  =  G  T 


where  we  have  pulled  the  dimensionless  time  1  out  of  the  temperature  in  the 

case  of  short  times  and  for  long  times  set  T  =  1.  For  the  anistropic  part  it  will 

3 

be  convenient  to  incorporate  the  factor  of  h  in  (23)  and  the  factor  of  rW2  in  (21) 
into  the  stress  optic  coefficients  of  equations  (17)  and  (18),  so  (23)  becomes, 

(25)  Ajfc  =  Af  -  AB.  t 

1,  2  o  1,  2 

The  anisotropic  part  is  characterized  by  the  two  constants 

(26) 

(27) 

Thus  the  whole  problem  is  characterized  by  three  material  constants. 

HI  DISCUSSION 

We  will  make  a  computer  parameter  study  of  the  three  material  constants 
including  a  comparison  for  G  along  with  those  of  the  slab  beam  reported  in  the 
2nd  quarterly  report.  We  will  also  look  at  the  birefringence  alone  as  well  as  in 
combination  with  the  isotropic  effect. 


b«  = 

b|  =  <x„2  h  To  ( (1-v)  p]2 -vpn>  n 
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In  the  next  quarter  we  will  have  available  a  general  heat  equation  program 
from  Oak  Ridge  National  Labs  so  we  may  compare  the  accuracy  of  our  approximate 
solution  of  the  heat  equation  with  the  exact  solution.  Equally  we  hope  to  either 
acquire  or  write  a  thermoelastic  computer  program  to  solve  the  thermoelastic 
problem  in  simple  disc  windows  of  isotropic  or  cubic  material  with  general 
temperature  distributions. 
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3.  DISCUSSION 


Additional  absorption  measurements  on  high  purity  melt  grown  GaAs  samples 
prepared  by  horizontal  bridgeman  and  float- zone  techniques  have  convinced  us  that 
there  is  a  lower  limit  to  the  absorption  coefficient  at  10.6|im  for  such  samples  of 
about  0  848  . 008Um.  The  similarity  of  the  absorption  spectra  in  the  9  to  11  |im 
region  for  all  samples  measured  is  very  striking,  considering  the  very  different 
materials  processing  techniques  r  ^presented  in  our  selection.  Our  temperature 
dependent  measurements  have  not  yet  been  completed. 

The  absorption  coefficient  of  one  of  the  Bell  and  Howell  GaAs  samples  has 
been  decreased  appreciably  by  annealing  as  anticipated  by  Prof.  Copley,  but  the 
reduced  absorption  spectrum  did  not  fall  below  the  lower  limits  observed  in  other 
higher  quality  samples. 

We  are  undertaking  accurate  measurements  of  the  dielectric  constant  £  of  GaAs 
as  there  seems  to  be  some  ambiguity  in  available  values  for  this  parameter. 

Prof.  Hellwarth  has  stressed  the  importance  of  accurate  knowledge  for  £  for 
comparison  of  theoretical  expressions  with  experiment. 

The  acoustic  surface  wave  technique  for  measuring  absorption  coefficients  has 
been  developed  to  the  point  where  preliminary  experiments  to  determine  the  10.6  iim 
absorbance  of  quartz  have  been  performed.  Extraction  of  the  absorption  coefficient 
from  the  data  requires  numerical  analysis,  now  in  progress.  Quartz  is  used  in 
these  initial  experiments  because  the  acoustic  wave  technology  is  well  developed 
for  this  material.  During  this  period,  a  device  has  been  fabricated  to  launch  and 
detect  acoustic  waves  on  alkali  halide  samples.  Samples  of  low  0KC1  with  high 
quality  surfaces  have  been  ordered  from  Janos  Optical  Corp.  (Vermont)  to  ex¬ 
pedite  this  project.  We  plan  to  develop  an  in-house  capability  for  preparing  al¬ 
kali  halide  surfaces  of  sufficient  quality  for  the  acoustic  wave  technique. 

KC1  crystals  have  been  grown  from  the  ultra  pure  material  described  in 
previous  reports.  The  apparent  quality  of  the  crystals  (by  visual  inspection)  is 
improving  as  experience  with  the  apparatus  accumulates.  Insufficient  data  on 
absorption  coefficients  has  been  obtained  for  these  samples  to  report  a  meaning¬ 
ful  number. 

The  experimental  work  on  high  temperature  transport  properties  of  CdTe, 
including  self  diffusion  of  Cd  and  Te,  and  high  temperature  Hall  effect,  has  been 
completed.  While  the  results  are  being  analyzed,  steps  are  being  taken  to  obtain 
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material  for  studies  of  the  electrical  and  optical  properties  at  lower  temperatures. 
The  material  will  probably  be  purchased  from  II- VI  Inc.  (Pennsylvania). 

Theoretical  work  during  this  period  included  the  development  of  further  checks 
on  the  approximate  multiphonon  absorption  theories  using  exact  moment  calculations 
reported  previously.  Some  of  this  work  was  reported  at  the  NBS-ASTM  Symposium 
on  Damage  in  Laser  Materials  in  May. 

The  superposition  approximation  for  acentric  stress  distributions  reported 
previously  has  been  exploited  during  this  period  in  a  calculation  of  the  influence 
of  stress  induced  birefringence  on  thermal  slewing  of  acentric  beams.  The  theory 
is  strictly  applicable  only  to  mechanically  isotropic  substances,  but  we  believe 
that  it  will  lead  to  reasonable  quantitative  estimates  for  the  size  of  the  effect  in 
cubic  anisotropic  substances  such  as  the  alkali  halides.  Difficulties  with  the 
computer  driven  plotter  at  our  computation  facility  prevented  us  from  including 
the  results  of  the  theory  in  this  report. 
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